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Abstract 

We investigate the Coulomb blockade resonances and the phase of the trans- 
mission amplitude of a deformed ballistic quantum dot weakly coupled to 
leads. We show that preferred single-particle levels exist which stay close 
to the Fermi energy for a wide range of values of the gate voltage. These 
states give rise to sequences of Coulomb blockade resonances with correlated 
peak heights and transmission phases. The correlation of the peak heights be- 
comes stronger with increasing temperature. The phase of the transmission 
amplitude shows lapses by n between the resonances. Implications for recent 

experiments on ballistic quantum dots are discussed. 
PACS numbers: 73.23.Hk, 73.23.Ps, 73.40.Gk 
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I. INTRODUCTION 



Quantum dots have been intensively investigated both experimentally and theoretically 
[l| in recent years. In this paper, we present a theoretical study of the correlations of 
conductance peaks and of transmission phases that have been observed in recent experiments 
on quantum dots in the Coulomb blockade regime. 

Quantum dots are small islands of electrons that are only a few hundred nanometers 
in size and typically contain a few hundred electrons. The spectrum of a quantum dot is 
determined by the Coulomb interaction of the electrons and by the external electrostatic 
confining potential. The confining potential and hence the size and shape of a quantum dot 
can be controlled by external gates. This makes quantum dots an ideal tool for studying 
the properties of finite systems of interacting fermions. 

Experimentally, the spectra of quantum dots have been measured using optical (far- 
infrared) spectroscopy and/or transport experiments. In the latter case, the quantum dot is 
coupled via tunnel barriers to external leads. The conductance measured at a finite drain- 
source voltage reveals the excitation spectrum of the dot whereas the linear conductance 
yields the addition spectrum of the quantum dot. Both the excitation and the addition 
spectrum are dominated by the classical Coulomb blockade effect: Large conductance peaks 
are observed when the dot potential is tuned in such a way that the number of electrons 
on the dot can fluctuate without any cost in energy. These peaks are nearly periodic in the 
gate voltage on the dot. At consecutive peaks the number of electrons on the dot changes by 
one. At values of the gate voltage located between the positions of the conductance peaks, 
electron transport through the dot requires a large charging energy. Hence, the current 
between conductance peaks is strongly suppressed, the remaining current being mostly due 
to virtual tunneling processes (co-tunneling regime) 0. 

Typically, metallic quantum dots have such a large density of states that the Coulomb 
blockade oscillations can be described by classical theory which ignores the discreteness of 
the spectrum. The situation is different in semiconductor dots. Here, the mean single- 
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particle level spacing A can be larger than the temperature kT. The regime A ^> kT, F, 
where T is the strength of the coupling to the leads, is the resonant tunneling regime. In 
this regime, each conductance peak is mediated by a single quantum state of the dot. The 
peak height of the conductance resonance is a direct signature of the wave function of the 
resonant state. 

Some years ago, Jalabert, Stone, and Alhassid developed a statistical theory of the 
Coulomb blockade |§ in the resonant tunneling regime. In order to explain strong fluctu- 
ations of the peak heights of neighboring conductance resonances, they assumed that the 
eigenstates of a quantum dot can be described by random matrix theory. According to this 
theory no correlations are to be expected for the peak heights of neighboring peaks. More- 
over, the distribution of peak heights is predicted to be universal, and only determined by 
the fundamental symmetries of the system. While two recent experiments [|]|| appeared to 
have confirmed these predictions, a more careful examination of the experimental data casts 
some doubt on the validity of random matrix theory in describing the physics of quantum 
dots. Deviations from the predictions of random matrix theory include the appearance of 
4-5 correlated transmission peaks ||, non Wigner-Dyson related distribution of conductance 
peak spacing |5|J7]], and reduced sensitivity to magnetic flux j|,[5],[7]|| . Even stronger correla- 
tions have been reported in very recent experiments using a quantum dot embedded in an 
Aharonov-Bohm ring || where both the peak heights and the phase of the transmission am- 
plitude could be measured. Strong correlations within sequences of more than 10 resonances 
were found. 

In the experiments of Refs. ||,|5]|7]|| , kT was approximately (0.1 — 0.5)A. When the 
resonant tunneling limit A 3> kT is not fully met, also neighboring eigenstates around the 
Fermi energy (rather than only a single state) contribute to the conductance peak. Extending 
the random matrix theory description to the regime A > kT, correlations in the peak height 



were found ]I0| . Those temperature-induced correlations seem sufficiently strong to account 
for sequences of up to 5 correlated peaks as reported in Ref. |J. However, temperature alone 
can clearly be ruled out as the source of the much stronger correlations found in Ref. ||. A 



mechanism different from and acting in addition to temperature must be the origin of these 
correlations. 

In this paper we propose such a mechanism. Our model involves certain geometry- 
specific assumptions and is therefore restricted in generality and universality. Nevertheless 
it may still pertain to experiments on nearly integrable ballistic dots (cf. Section [TIT]). 
Our mechanism is a synthesis of two approaches developed earlier 0,0]. In Ref. [jnj] a 
scenario was described for peak correlations at vanishing temperature. It was argued that 
deformation of the confining potential of the dot generically gives rise to avoided crossings 
of the single-particle levels on the dot. As a result of such crossings one and the same 
eigenstate of the dot may dominate a sequence of neighboring conductance peaks |TT| and 



thereby cause correlations. Although several arguments in support of this mechanism were 
given a quantitative study of the resulting correlations has not been presented yet. A 



different line of thought was pursued in Ref. [0. That paper aimed at explaining the "phase 
lapse" observed in the experiment of Ref. |J (a different theoretical discussion of the phase 
lapse was presented in Ref. ]I3"|). It was shown \Y1\ that the phase-lapse behavior as well 



as strong finite temperature correlations of conductance peaks arise if the dot supports one 
particularly well conducting state. This state would dominate a sequence of conductance 
peaks provided it remained within an energy interval of order kT around the Fermi energy. 

Here we study peak-height correlations and the transmission phase in quantum dots 
taking a deformed harmonic oscillator as a specific example. Using this model, we show 
that deformation of the confining potential leads to peak-height correlations, in keeping 



with the arguments of Ref. pLTjl . We identify those eigenstates of the quantum dot that 
are most strongly coupled to the external leads and which, therefore, support the bulk 
of the current through the dot. Peak-height correlations are strongest for sequences of 
resonances mediated by such states. We also investigate the influence of temperature on the 



correlations. As expected from the studies in Ref. ||12|| , we find that temperature leads to a 
marked increase of the correlations. Combining the effect of deformation and of temperature, 
we obtain sequences of up to 30 conductance resonances with similar transmission phases 



and similar peak heights. We demonstrate that this same mechanism may also account for 
the phase-lapse behavior observed in Ref. 0. 

The paper is organized as follows: In the next Section we introduce a model for a 
deformed quantum dot and describe how correlations can arise. In Section 3 we calculate the 
conductance with the help of a master equation. In Section 4 we investigate the transmission 
phase. This phase can be measured in Aharonov-Bohm type experiments containing a 
quantum dot. The last Section gives a summary and a discussion of the limitations of our 
model. 



The confining potential of a quantum dot is often defined electrostatically in terms of a 
split gate. As depicted in Figure 1, this is an arrangement of electrodes on the surface of a 
heterostructure. When a negative bias is applied to the gates, the two-dimensional electron 
gas located some 100 nm or so beneath the electrodes will be depleted. The barriers through 
which electrons can tunnel between the leads and the dot are denoted by A and B. A voltage 
V g applied to the plunger gate P controls the chemical potential on the dot. A change of 
V g not only changes the number of electrons on the dot but also distorts the confining 
potential of the two-dimensional electron gas in a substantial way, causing a deformation of 
the quantum dot [[□]]. 

We consider the standard Hamiltonian H for a quantum dot coupled to leads, containing 
the Hamiltonians H L and H R of the left and right leads, respectively, the Hamiltonian of 
the isolated quantum dot H D , and the Hamiltonian H T for tunneling between the leads and 
the dot, 



II. THE MODEL 



H 



H L + H R + H D + H T , 
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RT = E (yjfcof 4 + h. c .) + 2 (o* 4 + h - c -) • 

n,k n,k 

Here, e k ,R and e n are the energies and a^' R and c n the annihilation operators for single- 
particle states in the leads and in the dot, respectively. For the Coulomb interaction on 
the quantum dot we use the constant interaction model with N = J2n n n the number of 
electrons on the quantum dot. The tunneling matrix elements V^ R ^ involve the overlap of 
wave functions in the leads and in the dot and are given below. 

We model the confining potential as an anisotropic harmonic oscillator potential. A har- 
monic potential has been used previously in studies of quantum dots |l4j and, at least for 
small dots, is believed to be a fair approximation to the true confining potential. Although 
we use a specific model, most of our conclusions apply to any sufficiently smooth confining 
potential for which the Hamiltonian is nearly integrable. Then, the transverse and longi- 
tudinal modes in the dot are nearly decoupled, cf. Eq. (0). This condition is met when 
the matrix elements of the perturbation violating integrability are smaller than the mean 
single-particle level spacing. However, even mild disorder or boundary roughness violating 
this condition will modify our picture considerably. 

The energy eigenvalues e n for the quantum numbers n = (n x , n y ) are given by 

e n = E(n x ,n y ) = hu x (n x + -) + %uj y {V g )(n y + -) - aV g + Eq. (2) 

To describe the deformation, we assume that the oscillator frequency u) y (V g ) = u x (l — j(V g — 
Vq)) in the transverse direction y depends linearly on the gate voltage V g while the frequency 
uj x in the longitudinal x-direction is held fixed. The parameter a relates the overall depth 
of the dot's potential to the gate voltage. The constants Eq and Vo determine the number 
of electrons on the dot at zero deformation. 

The dependence of the single-particle levels on the gate voltage is shown in Figure 
The shell structure of the isotropic harmonic oscillator (V g = Vq) is clearly visible. It 
survives for small values of V g but is eventually destroyed by deformation. Each shell q is 
characterized by non-negative integer quantum numbers q = n x + n y . In each shell, there 
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are levels depending weakly (strongly) on V g , characterized by large (small) values of n x 
and small (large) values of n y , respectively. These are referred to as "flat levels" and "steep 
levels", respectively. A small deviation from integrability will change the level crossings 
shown in Figure into avoided crossings. For nearly integrable systems, the wave functions 
retain their character across avoided crossings. Flat levels are particularly stable, their wave 
functions change little with deformation (or gate voltage) and remain self-similar even after 



several avoided crossings jll 



The matrix elements V L , V R for tunneling from the left and right lead to the quantum 
dot are given by |1| 



(3) 

J X = X B 



dx x ' v dx 

Here ip^ is the wave function with wave vector k in the left (right) lead , and & nx ,n y 
is the wave function in the dot. The integration extends in the y-direction and xb is 



arbitrary but must be located within the barrier ||15|| . We restrict ourselves to the case of 
a single transverse channel in each lead. The nodes of the wave functions of flat (steep) 
levels with large n x (n y ) are predominantly carried by the x-component (y-component, 
respectively). Thus, the wave functions of flat levels extend much further into the barrier 
region and have considerably larger matrix elements V L ^ R '(k, n x , n y ) than those of the steep 
levels. This important property is illustrated in Figure It has immediate consequences 
for the conductance at finite temperature: The very same single-particle state can dominate 



different Coulomb blockade resonances seen at different values of V g ||1 1|| . We now explain 
this feature qualitatively, postponing a detailed discussion to later sections. 

At low temperature (kT <C A) and for small bias voltage, the transmission through the 
dot can be qualitatively obtained from the mean-field approximation for the dot spectrum. 
In this approximation, each single-particle energy q is replaced by the effective value £j = 
£i + UJ2j^i( n j) PHTTf - According to Koopmans' theorem, is the energy needed to add 



an electron in state % to the dot, whereas the excitation energy at fixed electron number is 
given by the difference of the corresponding two effective energies. Because of the Coulomb 



interaction, there is a gap of magnitude U between the last occupied and the first empty 
effective single-particle level, while the other occupied (empty) levels below (above) the 
Fermi energy Ep are on average separated by the usual mean level spacing A. Avoided 
crossings of single-particle levels result in avoided crossings of the effective levels, in spite 



of this gap [IT! 



A Coulomb blockade resonance occurs, and the number of electrons on the dot changes by 
one, whenever an effective single-particle level crosses the Fermi energy Ep of the reservoirs. 
(We assume Ep to be independent of the gate voltage in the following). For U 3> A, the 
distance between adjacent resonances is 5V g = U/a. Without level crossings, different 
resonances correspond to different single-particle levels. In the presence of level crossings, 
the situation changes. This is shown in Figure |] which displays the gap between the filled 
levels below and the empty ones above Ep. Resonances occur at gate voltages V\,Vi and 
V3. Suppose a flat level F (dashed) crosses Ep at V\. If there is an avoided crossing of 
F with a steep level Si from a higher shell between V\ and V2, level F is pushed above 
Ep while level S\ is immersed into the Fermi sea. At V2 the flat level F crosses Ep again, 
causing another resonance to occur. The mechanism works again between V2 and V3 where 
another steep level S% intersects with F. Since flat levels keep their wave functions after 
avoided crossings, the resonances at V\,V2 and V3 all carry the same single-particle wave 
function. This mechanism gives rise to strong correlations of the properties (peak height 
and transmission phase) of several resonances. However, it leads to strong correlations only 
if there is one and only one crossing of the flat level F with a steep level within subsequent 
intervals 8V g . If - on average - there is less than one (more than one) crossing, the level F 
will eventually be pulled down below (up above) the Fermi energy and will become irrelevant 
for the behavior of Coulomb blockade resonances. 

So far, we have focussed attention on the resonant tunneling regime T, kT <C A where 
only the level at the Fermi energy determines the properties of the resonance. Essential 
modifications arise for finite temperatures kT ~ A. Here, also levels at a distance ~ kT 
from Ep contribute to the resonance. Since the flat levels are coupled to the leads much 
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more strongly than the steep ones, the presence of a flat level at a distance ~ kT from 
Ep suffices for it to dominate the resonance. Long sequences of correlated resonances may 
occur if repeated avoided crossings cause a flat level F to stay sufficiently close to Ep over a 
sufficiently long range of V g values. This is the picture we investigate quantitatively for the 
case of an anharmonic oscillator in the sequel. The picture suggests that the correlations of 
consecutive Coulomb blockade resonances will increase with temperature. 

The number of intersection points of a flat level (n x ^ , n y = 0) with steep levels from 
higher shells can be determined from Eq. (fj). The total number of crossings of the flat level 
occurring in the range V g — Vq is given by 

C 2 1- 7 (V,-Vo)' [) 
The number of steep levels from higher shells increases with deformation and causes N c to 
increase, too, until it diverges for the unphysical situation of extreme deformation uo y — > 0. 
In order to have one crossing within an interval 5V g = Ufa we require (dN c )/(dV g ) = a/U. 
This condition yields the value of V g where maximal correlations between Coulomb blockade 
resonances should occur. We estimate the number AN of resonances for which the distance 
between the flat level (n x , 0) and the Fermi energy is less than one level spacing and obtain 



AN ~ 2, 



with a deformation u) v /uj x = J(jU)/(2a) n x . 

III. COULOMB BLOCKADE RESONANCES 

In this section we calculate the conductance G of a quantum dot with the single-particle 
spectrum (fj) in linear response. We use the master equation |18[]. In the Coulomb blockade 
regime, maxima in G occur for values of V g where the configurations with iV and N + 1 
electrons on the dot are degenerate. The resulting sharp conductance peaks are almost 
equally spaced PHTB|. For kT A each peak is due to a single level n of the dot, and 
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the peak height is given by G ~ T^T^/(T^ + T%) @ with the tunneling rates ~ 
J2k \Vkn\ 2 $( e n ~ e t^)- We note that the peak heights are highly sensitive to the wave 
functions in the dot. Strongly coupled levels give higher peaks than the weakly coupled 
ones. For finite temperature several single-particle states contribute to a resonance. 

Under the assumption of sequential tunneling, transport through the dot at finite tem- 



perature is described by the master equation ||18|| . In the regime kT 3> T this equation 
determines the occupation probabilities P v of the single-particle levels of the dot under the 
influence of the interaction U and of the coupling to the leads. It is given by 

d 

-^Pu = p i£in{n -^v)- P v Y out {y -> fj) . (6) 

Here \x and v label Fock states, i.e., Slater determinants defined in terms of the occupation 
numbers of all single-particle states in the dot. The symbols T in and T out stand for the rates 
of the tunneling processes into and out of the dot. These processes change the number of 
electrons on the dot by one, and the associated Fock states from /i to u, and vice versa. The 
rates contain not only the coupling of the specific single-particle states of the dot to the 
leads but also take into account a possible suppression of tunneling by the occupation of the 
states in the leads. By expanding Eq. (O) around the equilibrium probability distribution 



Beenakker [18] obtained the conductance G for small bias voltage, 



p 2 oo -pL-pR 

G=7^EE flXT^ P ^ t 1 - F ^\N)\ f(e n + U-N). (7) 

Here, / is the Fermi function and P^ q is the probability to find N electrons on the dot, 

= tr N exp(-l3H D ) ti N exp(-pH D ) 
N trexp(-/3H D ) E7vtr JV exp(-^) ' 1 ' 

The inverse temperature is denoted by /3, and tr^ denotes the trace over the Fock states 
with iV electrons on the dot. The canonical occupation number of level n when there are N 
electrons on the dot is given by 

F (e i \N) = (n i )N- tlNexp{ _ pHD) ■ (9) 
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Eqn. (0) has also been derived Jl9| using a Landauer-Biittiker type approach generalized 
to include the interaction of electrons on the dot. 

Numerically, it turns out to be sufficient to calculate F eq (e n \N) for a window of levels 
around the Fermi energy, and to take the occupation numbers equal to 1 below and 
above this window. The following results are obtained using a window of 16 levels that are 
populated with 8 electrons. We checked that our results are insensitive to changes of the 
window size. 

Using Eq. (f|) and the relation T^ RS> ~ J2k iVfcn | 2 ^( e n — e t )> we ^ n( ^ that the tunneling 
rates T n are proportional to the square modulus of the harmonic oscillator wave functions 
at the position of the barriers f— (xb, Vb), 

\2 



2^+ n v)n x \n y \ ' K } 

where H n are the Hermite polynomials and where we have used dimensionless variables 



ix = ymLu x /h xb and £ y = ymujy/h ys- Neglecting deformation we can relate the barrier 
height Vb to the position of the barrier via Vb/TiuJx — (£| + £y)/2- In the sequel we choose a 
fixed value Vs/huj = 25 for the barrier height. [] We consider two different geometries of the 
leads connecting the quantum dot with external reservoirs: (i) The two leads are located 
exactly opposite to each other, so that £ x = ±v^50 and £ y = 0. (ii) The leads are arranged 



at an angle of 90 degrees, with the barriers at £ x — £ y — — y50/2 and £ x = —£ y = y50/2. 
The latter geometry has been used, for instance, in the experiment of Ref. 0. 

The tunneling rates T n for both geometries and for several states n are presented in 



Compared with typical experimental parameters this value appears to be too small. We use this 
value in order to avoid an unphysically large increase of V with increasing quantum number n x . 
Such a strong increase is characteristic of the harmonically shaped barrier. A less steep increase 
would be obtained for steeper tunneling barriers. Such barriers appear to give a more realistic 
description of the depletion zone of the electron gas near the gates. For reasons of consistency with 
the harmonic oscillator model used throughout this paper, we decided to use harmonic barriers. 
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Table |. We note that for geometry (i) the flat levels are coupled much more strongly to 
the leads than the steep ones (columns 1 and 2). Comparing the rates for the flat levels of 
different shells we find a considerable increase with increasing n x (columns 4 and 1). For 
geometry (ii) the states within one shell with equal quantum numbers n x and n y are most 
strongly coupled to the leads. However, the difference between strongly and weakly coupled 
states is not as pronounced as in geometry (i). 

We conclude that independently of the precise shape of the barrier, in a geometry of type 
(i) which is realized in Figure [l| flat levels are more strongly coupled to the leads than steep 
levels. This is because the wave functions of flat levels have cigar-like shapes extending 
closer to the leads (cf. Figure |3|). Therefore, flat levels carry the bulk of the current. 
The difference between well coupled and poorly coupled dot states is less pronounced when 
the leads are arranged at an angle (geometry (ii)). In this case, the mechanism for peak 
correlations described in Sec. |H| leads us to expect a reduction in the length of sequences of 
correlated peaks. 

We present results for a = 1, 7 = 0.005, E Q = —11 and V = 90. At zero deformation, 
the mean level spacing A is related to the harmonic oscillator frequency u x in x-direction 
by 

. _ -£f _ hu x N sh ^ 2hu x 

~ N el ~ l(N sh + l)(N sh + 2) ~ N sh ' 1 j 

Here N s h is the number of the last filled shell and N e i = (N s h + l)(N s h + 2)/2 is the total 
number of electrons on the dot. We take N s h = 14 so that there are about 100 electrons 
on the dot. We assume A = 0.03C7 which roughly corresponds to the situation of the 
experiments of Refs. PH . 

In order to monitor the influence of a flat level F on the conductance, we define the 
distance d of F from the Fermi energy as the number of levels between F and the Fermi 
energy including F itself and count positively (negatively) for states above (below) the Fermi 
energy. Thus, d = 1 indicates that F is the first unoccupied level. Figure [5] shows d vs. gate 
voltage for the flat level n x = 14, n y = 0. Adjacent points correspond to adjacent Coulomb 
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blockade resonances. For 130 <V g < 160, F stays in the vicinity of the Fermi energy. Hence 
there is - on average - one crossing with a steep level in the interval 5V g . For V g < 130 
(V g > 160), the number of avoided crossings in the interval SV g is less than one (bigger than 
one), and F moves towards the (away from the) Fermi energy, respectively. The jumps at 
V g = 156 and V g = 180 are due to multiple crossings of levels which occur because of the 
integrability of our model. 

For the same choice of parameters as in Figure [| and for geometry (i) , Figure |6| shows 
the conductance vs. gate voltage for two temperatures, (a) kT = 0.2A and (b) kT = 0.4A. 
About 100 Coulomb blockade resonances occur in the interval 100 < V g < 200. In both plots 
strong peaks with similar peak heights appear whenever the flat level F is close to the Fermi 
energy, especially at the higher temperature (case (b) of Figure |6[). In the regions V g < 130 
and V g > 160 where the conductance is not dominated by the flat level F, G is much smaller 
than in the interval 130 <V g < 160. On the scale of Fig. ||, some of the conductance peaks 
are not even visible. 

Figure [7] shows the conductance vs. gate voltage for the same parameters as in case (b) 
of Figure ^ but for geometry (ii). Now, the resonances are dominated by steeper levels from 
higher shells which are coupled more strongly to the leads than the flat level. This is why 
the peak heights are bigger on average than in Figure |6|, why they show stronger variation, 
and why they increase systematically with increasing gate voltage. 

IV. PHASE 

We now turn to the behavior of the phase of the transmission amplitude through the 
quantum dot. This phase has recently been measured in a set of experiments using an 
Aharonov-Bohm (AB) interferometer with a quantum dot embedded in one of its arms. We 
consider the simplest case where the AB interferometer is coupled to only one channel in 
each connecting lead. The transmission coefficient T through the AB device is then given 
by 
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T « To + 2Re jt*e- 2 ™*/* 1 d£ ^-J^ t QjD (,E7) j . (12) 

Here, % = | to | 2 is a flux- and energy-independent term given by the square of the amplitude 
for transmission through the empty arm of the AB interferometer, while tQD is the amplitude 
through the arm containing the dot. Since the quantum dot is weakly coupled to the arm, 
we have |£q_d| <^ l^ol- We have explicitly displayed the dependence on the magnetic flux <E> 
through the AB device and neglected higher harmonics. The symbol $o = h/e denotes the 
elementary flux quantum. 



The master equation used in Section |J| deals with occupation probabilities and is, 



therefore, not able to yield the phase of the transmission amplitude tq^. We have used 
another approach. We have expressed tqu in terms of the retarded Green function G ret of 
the dot, 



t QD (E) = £ V l L (E)G^(E)Vf*(E) . (13) 

The finite-temperature Green function G ret must be calculated in the presence of the inter- 
action U and the tunneling. A derivation of G ret starting from the equations of motion is 
given in Appendix [A]. Assuming that the total number of electrons on the dot is a constant 
of motion, we obtain 



oo 



N=0 



1 - (Wi)jV , ("i)jV 



E-(ei-fi+ UN) + iY t /2 E-{ei-n + U(N - 1)) + ir</2 

Here Tj = Vf + Tf . The probability P^ q that there are iV electrons on the dot and the 
canonical occupation number (n)jy are given in Eqs. (§) and respectively. 

Within our approximations the Green function G ret is diagonal. This fact implies that 
real ( particle-hole) excitations of the dot caused by tunneling transitions are not taken into 
account. This is justified in the regime of elastic cotunneling kT < \fU~K where inelastic 
cotunneling processes do not contribute significantly to the transmission 0. Equation (|I4]) 
is a good approximation to the exact retarded Green function between Coulomb blockade 
resonances where fluctuations in the occupation number of the dot are strongly suppressed. 
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• (14) 



Moreover, even at resonance where G ret reduces to a single Breit-Wigner term, Eq. (|14]) is 
expected |^IJ to apply provided there are no degeneracies and we work well above the Kondo 
temperature pifl . The success of Eq. (|TJ|) in these limiting cases suggests that well above 
the Kondo temperature, Eq. ([14]) is a good approximation to the exact Green function for 
all energies. 

Combining Eqs. (0) and (|i~2l) we obtain 



T = T + Ret* Q t QD e- 2 ^ , (15) 



where 



l QD 



dE {-%) tQD{E) 







Z7r< i N=0 



(^))^ (2) (£-.{Ki-<* + »-UN) + \ 
\ 2m 2 



(hi) ^ (2) (^-m -ei + fi-U(N-l)) + ^ 



(16) 



with the trigamma function ijj^ 2 \ 

In Figure || we show the phase <fi of the transmission amplitude versus gate voltage. As 
in the calculation of the conductance in Section |TTT| , the canonical occupation numbers are 
obtained by distributing 8 electrons over a window containing 16 levels. We take Tj = T = 
0.002f/ = A/15. The solid lines at the bottom of the plots show the conductance peaks and 
help to identify the resonance positions. In the left part of Figure || the flat level n x = 14, 
n y = is close to the Fermi energy. Here we find a strikingly similar behavior of the phase 
at all resonances. This behavior is found not only within the V g interval shown but for the 
entire interval 130 < V g < 160 comprising 30 resonances. The phase regularly increases by 
7r at resonance and displays a sharp lapse by 7r between adjacent resonances. As observed in 
Refs. (12|,[13j the increase at resonance occurs on the scale kT (we assume kT > T) and the 
phase lapse between resonances on the scale T. The temperature dependence of the phase 
is shown in Figure ^. In the right part of Figure |8] we show the transmission phase for the 
case where the distance between the flat level and the Fermi energy is large compared to 
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kT and increases with V g (cf. Figure |5|). The phase behaves less regularly, with an increase 
by 7i or less at and an immediate phase lapse near the resonances. Between resonances the 
phase remains virtually constant. 

To interpret our results, we consider first the phase <ft at resonance. The identical behavior 
of 4> at all resonances in Figure |8|(a) reflects the fact that at each resonance, the transmission 
through the dot is dominated by the strongly coupled level F. This is the same mechanism 
as in the sequence of strong conductance peaks shown in Figure @(a). The more erratic 
phase behavior seen in Figure §(b) is the result of the interplay of various levels of the 
dot. The regular behavior of the phase lapse between adjacent resonances is also due to the 
dominance of the flat level. At finite temperature the flat level F has a finite probability 
of being either occupied or empty and, thus, may contribute to both an electron-like and 
a hole-like cotunneling process. The contribution of both processes to the transmission 
amplitude through the dot is 



L\rR* 



F v F 



1 - (n F ) N (n F ) N 



(17) 



_E - (e F - fi + UN) + iT F /2 E - (e F - (i + U(N - 1)) + iT F /2_ 

where the first (second) term represents the electron-like (hole-like) contribution, respec- 
tively. As the gate voltage V g = nja scans the N th valley (i.e. varies from (ejy+U-(N— l))/a 
to (ejv +1 + U ■ N)/a), the sign of Re t F reverses, leading to a phase lapse. The same conclu- 



sion has previously been reached in Ref. [12|; an interpretation in terms of scattering theory 



has been given in Ref. [13|. If F is far away from the Fermi energy (on the scale of kT) 
either the particle-like or the hole-like process will dominate, and the phase lapse moves 
from the valley towards the resonance, as depicted in Figure |8|(b). 

We emphasize that the phase lapse between resonances is a genuine interaction effect. 
Indeed, the interaction U is needed to keep the flat level close to the Fermi energy for a long 
sequence of resonances. For non-interacting particles (U — > 0) the transmission amplitude at 
different resonances would be dominated by different single-particle levels. In the same limit, 
the cotunneling amplitude flTTD would reduce to a single, temperature-independent term. 
The phases of the transmission amplitude in consecutive valleys would not be correlated, 
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and there would be no systematic phase lapse between resonances. We also note that the 
systematic phase lapse occurs only at finite temperature. At zero temperature a flat level 
could only contribute to either particle-like or hole-like cotunneling. 



V. SUMMARY. THE QUESTION OF NON-UNIVERSALITY 

Since the first measurements on quantum dots in Aharonov-Bohm interference devices 
were reported, many aspects of phase-coherent transport through quantum dots have been 
understood theoretically. However, one of the most striking features, the strong correlations 
of the transmission phases in sequences of many resonances, has long withstood a satisfactory 



theoretical explanation. Earlier attempts [0,0 to solve the problem could account for short 



sequences but not for the sequences of more than 10 resonances found experimentally. 
In this paper we have demonstrated the viability of a mechanism, based on a synthesis 



of the ideas proposed in Refs. ||n| , [Llq| , that gives rise to long sequences of correlated peak 
heights and transmission phases. We have used several approximations, the most central one 
being that the confining potential defining the dot is "almost" integrable. More precisely, 
both the deviation from integrability and the disorder must constitute a perturbation which 
is small on the scale of the mean single-particle level spacing. Our model shares many 
features of the ballistic quantum dots employed in some experiments PU] and may account 
semiquantitatively for some of the experimental observations. Nevertheless, our analysis 
naturally falls short of providing a complete and universal framework which could account for 
the combined effect of disorder and interaction on correlations in transmission experiments. 

We have used the following specific conditions and assumptions: (i) Among the eigen- 
states of the quantum dot, some must be coupled more strongly to the leads than others. 
This assumption is met by a model which is nearly integrable, as is the case for our parabolic 
confining potential. This potential renders (n x ,n y ) good quantum numbers for all values 
of V g . Of all levels in a shell, the level with n y = is most strongly coupled to the leads, 
(ii) Changes in the gate voltage induce deformations of the dot boundary such that the po- 
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tential is deformed in the transverse ^/-direction. This assumption guarantees that in each 
shell, the level with n y = is flat, i.e., stays close to the Fermi energy over a wide range 
of Vg-values. (iii) On average there is one crossing of the flat level with one other level per 
unit interval. This interval is defined by the change in V g needed to add an extra electron 
to the dot. (iv) The temperature is sufficiently high to produce sufficiently long sequences 
of correlated resonances. The minimum temperature required by this condition depends 
both on the distance of the most strongly coupled level from the Fermi energy, and on the 
relative strength of the coupling of that level to the leads. With increasing temperature the 
correlations become more robust. 

Strong boundary deformations leading to strongly chaotic classical motion within the 
dot, or strong disorder in the dot are likely to destroy the correlations altogether since they 
generically do not allow for the existence of eigenstates that are particularly well coupled 
to external leads. In this sense, the correlations proposed in the present paper are non- 
universal in origin. This conclusion agrees with the observation of weak peak correlations 
in the strongly deformed dots of Ref. || as compared with the strong correlations found in 
the experiments of Ref. || . 

Our ideas may be checked experimentally on dots that are not embedded in an AB device 
but are coupled directly to leads. This setup does not allow for tests of phase correlations but 
provides a convenient setup for measuring conductance peak correlations. The conductance 
of a dot with a regular (rectangular) lithographic shape has been measured by Simmel et al. 
0. These authors did indeed find a sequence of more than 10 strong peaks with very similar 
peak heights. In the same sweep they also observe envelopes of smaller peaks very similar 
to our results in Figure || To test our picture further, it would be illuminating to perform 
a similar two-terminal conductance experiment with leads attached at two sides of the dot 
which form an angle of 90 degrees. Here, a reduction of the correlations is to be expected. 
It would also be interesting to compare two setups, one with a plunger gate and the other 
with a backgate configuration. In the latter case the potential deformation is reduced. This 
should suppress our correlation mechanism. 
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APPENDIX A: RETARDED GREEN FUNCTION 

We derive Eq. (fH]). The retarded Green function is defined as 

Glf = - iO(t) ([c k (t), 4(0)]+) = F kl + G kl where (Al) 
Gki = - i9(t)(c k (t)cj(0)) F kl = -i9(t)(cj(0)c k (t)}. (A2) 

The brackets denote the thermal average, (...) = tr(... exp(— /3if))/£r(exp(— /3H)). With 
Pn = trw exp(~(3H)/trexp(—(5H), we write the trace as a sum of terms with a fixed 
number N of electrons on the dot, 

oo oo 

G kl = -i9{t) £ P N (c k (t)4(0)) N = £ Gg° • (A3) 

N=Q N=0 

In the "equation of motion" method the Green function is differentiated with respect to time. 
Since the time evolution of an operator is given by the commutator with the Hamiltonian, 
a system of differential equations containing higher-order Green functions is generated. A 
closed system is obtained if these Green functions can be approximately uncoupled. The 
solution is obtained by Fourier transformation. Specifically, 

J^f ) = -i5(t)( Ci (0)c}(0)) N - 9(t)([ Ci ,H}(t)c](0)) N (A4) 

= -i8(t)(ci(p)c](p)) N - 9(t) ((3 - fj,)(ci(t)4(0)) N + U(N(t) Cl (t)c](0)) N ) 

-d(t) E (v^Mt) L c}(0)) N + V k T(a k (t) R c}(0)) N ) . (A5) 
k 

Since the interaction contains two creation and two annihilation operators a two particle 
Green function appears in the second line. The last two terms stem from the coupling to the 
leads, and by another equation of motion can be expressed in terms of the Green function 
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for the dot. Assuming that the states in the lead and in the dot are uncorrelated at t = 
we obtain 

(a k (t)c}(0)) N = J dlG^\t - t)6(t) exp(-ie£t) . (A6) 

i 

Fourier transformation of Eq. ( |A5| ) yields 

cuGlf } = (1 - m)^ + (ei - fi)G\P + UG i:j + 

V ki V kl , ^fci "feZ \ ^(JV) / A 

w-et + iS + ^WTiS) j (A?) 

where G~ij = —i J dt9(t)(N(t)ci(t)cj)^ expiout. The equations of motion are closed by assum- 
ing the total occupation A" to be constant. Then the number operator can be taken out of 
the expectation value. The assumption is justified in the valleys between resonances whereas 
at each resonance, A" fluctuates. For isolated resonances (level width <C level spacing), we 
have Im £ fc V£*V$/(u - t L k + iS) = -i5 u Tf/2. This yields 

r (N) M 1 -nj) N 

ij u - (ei - ii + UN) + iTf/2 + iVf/2 ' 1 } 

For Fij we proceed analogously and eventually obtain Eq. (|14]). 
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TABLES 



n = (n x ,n y ) 


(14,0) 


(0,14) 


(6,8) 


(12,0) 


(6,6) 


T n for geometry (i) 


2.68 • 1(T 6 


4.04 • 10" 23 


5.34 • 10~ 14 


9.09 • 10~ 8 


6.11 • 10~ 14 


T n for geometry (ii) 


6.30 • 10~ 12 


6.30 • 10~ 12 


5.66 • 10~ 7 


2.88 • 10~ 12 


2.41 • 10~ 8 



TABLE I. Tunneling rates for several levels n calculated from Eq. (|lG|). In geometry (i) the 
leads are opposite to each other with £ x = \/50 and £ y = 0. In geometry (ii), the leads are arranged 



at an angle, with £ x = ^ y = ^/50/2. 
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FIGURES 
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FIG. 1. A quantum dot (schematic). Electrons enter or leave the dot through the tunneling 
barriers A and B. The plunger gate P controls the number of electrons on the dot. Equipotential 
lines of the confining potential are shown for the isotropic (dashed) and deformed (dotted) case. 

^ shell q 




V Q Gate Voltage \ g 



FIG. 2. Dependence of the single-particle energies of the dot on the gate voltage V g . 
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FIG. 3. The thin solid line shows a cross section of the potential in longitudinal direction, the 
two barriers lying at opposite ends. The overlap of the dot wave functions (probability shown as 
dashed lines) and of the lead wave function (probability shown as a solid line on the left) increases 
strongly with the quantum number n x in x-direction. 
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E 




FIG. 4. Avoided crossings with the steep levels S± and S2 cause the flat level F (dashed) 
to stay close to the Fermi energy Ef- Resonances dominated by this level occur at gate voltages 
Vi, V2, V3. There is a gap of magnitude U between the last occupied and the first empty level. 
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FIG. 5. Distance d (in number of levels) of the flat level [n x = 14, n y = 0) from the Fermi 
energy Ep versus gate voltage V g . 
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FIG. 6. Conductance G vs. gate voltage for (a) kT = A/5 (left) and (b) kT = 2 A/5 (right). 
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FIG. 7. Conductance G vs. gate voltage V g for the same parameters as in Figure || (b) but 
with the leads arranged at an angle of 90 degrees. 
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FIG. 8. Phase (j) of the transmission amplitude versus gate voltage V g at kT = A/5 in two 
different intervals. The solid lines at the bottom of the plots display the conductance peaks. In 
the case shown in the left (right) part, the fiat level n x = 14, n y = is at or near (far removed 
from) the Fermi level, respectively. In the case of the right part, the flat level influences the phase 
as a background only. 
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FIG. 9. Transmission phase for kT = 0.2A (open circles) and kT = 0.4A (filled circle). The 
increase by tt at the resonance takes place on the scale kT, the phase lapse between resonances, on 
the scale T. 
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